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ABSTRACT 

We propose to use multiple-imaged gravitational lenses to set limits on gravity theo- 
ries without dark matter, specificly TeVeS (Bekenstein 2004), a theory which is con- 
sistent with fundamental relativistic principles and the phenomenology of MOdified 
Newtonian Dynamics (MOND) theory. After setting the framework for lensing and 
cosmology, we derive analytically the deflection angle for the point lens and the Hcrn- 
quist galaxy profile, and study their patterns in convergence, shear and amplification. 
Applying our analytical lensing models we fit galaxy-quasar lenses in the CASTLES 
sample. We do this with three methods, fitting the observed Einstein ring sizes, the 
image positions, or the flux ratios. In all cases we consistently find that stars in galaxies 
in MOND/TeVeS provide adequate lensing. Bekenstein's toy [i function provides more 
efficient lensing than the standard MOND ft function. But for a handful of lenses a 
good fit would require a lens mass orders of magnitude larger/smaller than the stellar 
mass derived from luminosity unless the modification function fi and modification scale 
do for the universal gravity were allowed to be very different from what spiral galaxy 
rotation curves normally imply. We discuss the limitation of present data and sum- 
marize constraints on the MOND /i function. We also show that the simplest TeVeS 
"minimal-matter" cosmology, a baryonic universe with a cosmological constant, can 
fit the distance-redshift relation from the supernova data, but underpredicts the sound 
horizon size at the last scattering. We conclude that lensing is a promising approach 
to differentiate laws of gravity (see also astro-ph/05 12425). 

Key words: gravitational lensing — cosmology, gravity 



1 INTRODUCTION 

The standard paradigm of Einsteinian gravity with dark 
matter and dark energy has proven amazingly successful at 
describing the Universe, in particular the data from the Cos- 
mic Microwave Background (e.g. Spergel et al 2003), galaxy 
redshift surveys (Percival et al 2002, Tegmark et al 2004), 
Type la supernovae (Reiss et al 1998, Perlmutter et al 1999), 
and weak gravitational lensing (see e.g. van Waerbeke & 
Mellier 2003, Refregier 2003 for reviews) to high accuracy 
with a small number of parameters. However, it is worth 
exploring alternative models of gravity to assess the unique- 
ness of the model and to open up new ways to unify gravity 
and the standard model. Indeed, the detection of deviations 
from the Einstein-Hilbert action could signal new physics, 
and such deviations are expected, for example, in models 
such as the M-theory inspired braneworld. 

The central role of both dark matter and dark energy in 
the cosmological model has also led some workers to ques- 



tion the standard paradigm. Given this it seems sensible 
to develop methods to test the basic assumptions of this 
paradigm, if only to put them on a firmer basis. One of the 
most direct ways to probe gravity over large scales in the 
Universe is via the effect of gravitational lensing. 

In this paper we shall explore gravitational lensing 
in the recently developed relativistic version of Modified 
Newtonian Dynamics (MOND), the Tensor- Vector-Scalar 
(TeVeS) theory, developed by Bekenstein (2004). The non- 
relativistic version of MOND was originally proposed by Mil- 
grom (1983) as an alternative to the dark matter paradigm. 
Milgrom (1983) suggested that galaxy rotation curves V(r) 
could be explained by modifying gravity: 
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where the interpolation function fi(g/ao) is the effective "di- 
electric constant" , which itself has the above asymptotic de- 
pendence on the gravitational field strength, g. Thus the 
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gravitational field strength g becomes significantly stronger 
than Newtonian gravity GM/r 2 in the weak MOND regime, 
when g < ao. 

From a theoretical point of view, MOND has not just 

one extra free parameter, ao, which is tuned to explain 

Tully-Fisher relations, but also a free interpolation function, 

(i(g/ao). With a standard choice n(x) = , x , one can 

Vl+i 2 

achieve very good fits to contemporary kinematic data of a 
wide variety of high and low surface brightness spiral and 
elliptical galaxies; even the fine details of velocity curves 
are reproduced without fine tuning of the baryonic model 
(Sanders & McGaugh, 2002; Milgrom & Sanders, 2003). 
However, whether MOND qualifies as a gravitational the- 
ory depends on its prediction of fundamental properties of 
gravity, e.g., the bending of the light. 

In a more rigorous treatment of MOND (Bekenstein 
& Milgrom, 1984), gravity is the gradient of a conserved 
potential, $(r). It satisfies a modified Poisson's equation, 
and trajectories of massive particles are governed by the 
equation of motion as follows, 

^=g = -V$(r), V. [g^g/ao)] = -AirGp, (2) 

where the right-hand equation is the density of all baryonic 
particles, g = |g| is the magnitude of the gravitational field. 

These versions of MOND, however, all suffer from be- 
ing non-relativistic. The main problem is that the theory is 
not generally covariant, as the physics still depends on the 
measured local acceleration. In the absence of a relativistic 
version of MOND, the paradigm could not be used to build 
cosmological models and could not provide robust predic- 
tions for the expanding Universe, the Cosmic Microwave 
Background, the evolution of perturbations, and gravita- 
tional lensing. 

Recently a fully relativistic, generally covariant version 
of MOND has been proposed, TeVeS (Bekenstein 2004), 
which passes standard local and cosmological tests used to 
check General Relativity. In this relativistic version, the con- 
formal freedom of a general relativistic model is used, along 
with a new scalar field, one vector field and a conformally 
coupled metric tensor, to preserve general covariance. Fol- 
lowing Bekenstein's (2004) paper, a number of other works 
have appeared studying the cosmological model (Hao & 
Akhoury, 2005) and large-scale structure of the Universe 
(Skordis et al., 2005) in the relativistic TeVeS theory. There 
are also attempts to refine the Lagrangian for the scalar 
field in TeVeS (hence the MOND interpolation /i-function) 
to improve fits to galaxy rotation curves (Famaey & Binney 
2005) and avoid unphysical dilation effects in two-body sys- 
tems (Zhao & Famaey 2005). A further extention of TeVeS, 
dubbed BSTV, has also been recently developed by Sanders 
(2005), which adds more flexibilities into the theory by mak- 
ing one of the non-dynamical scalar field in TeVeS a dynam- 
ical field. It remains to be seen what the range of predictions 
of these theories are in cosmology and solar system dynam- 
ics. 

Our goal here is to develop and test TeVeS by making 
concrete predictions of lensing in the theory, and comparing 
these predictions to data from strong lensing. We present 
the general framework for lensing in TeVeS. As a first ap- 
plication to data, we apply models of simple axisymmetric 
lenses: the point-mass lens and the Hernquist profile. While 



the point-mass is a very poor model for galaxies in a dark 
matter theory, it is a reasonable model for lensing in TeVeS, 
where the baryons are concentrated in a central region, and 
the Einstein ring in most strong lenses encloses most of the 
baryons of the galaxy. As an extension to this, we also con- 
sider the Hernquist profile, which allows us to develop our 
analysis to include extended galaxies. 

Our lensing formulation bears similarities with that of 
Qin et al. (1995), who developed a heuristic formulation to 
compute lensing in MOND in the weak field regime. Here we 
develop a more rigorous approach based on the relativistic 
TeVeS. In the weak-field regime of Bekenstein's model we 
find a bending angle a factor of two greater than that found 
by Qin et al. (1995), consistent with GR. 

As reviewed by Bekenstein (2004) and Sanders (2005), 
there are probably many ways of generalising MOND to a 
relativistic theory, with TeVeS (or BSTV) being the most 
successful one so far. Although we derive our results within 
the framework of the TeVeS theory, we explore a regime 
where dynamics behaves like in MOND. Therefore our re- 
sults are likely relavent for other relativistic versions of 
MOND as well. 

Nevertheless fundamental differences exist between 
MOND and TeVeS, not only at the conceptual level, but 
also at the phenomenology level, particularly when dealing 
with a non-spherical potential with tide or external field. 
For example, the Roche lobe of a satellite in MOND is more 
squashed in MOND than naive extrapolations of the New- 
tonian result (Zhao 2005), and the shape depends on the 
MOND function fi (Zhao & Tian 2005) or the TeVeS ^ for 
the scalar field (Zhao & Famaey 2005). Major development 
is also made in numerical tools to solve MONDian Pois- 
son's equations for studying systems of a general geometry 
in detail (Ciotti et al. 2005). For these reasons, we prefer 
to follow the TeVeS formulation of lensing and cosmology 
instead of previous phenomenology-based formulations (Qin 
et al. 1995, Mortlock & Turner 2001). 

The paper is set out as follows. After a brief presen- 
tation of the basic equations of Bekenstein's TeVeS theory 
in §2. We give in §3 the central equations for gravitational 
lensing, and show how to calculate gravitational potential 
around a galaxy. An eager reader could go directly to §4, 
where we model distances and ages in a homogeneous and 
isotropic, expanding universe, and determine the parameters 
for the minimal-matter cosmology by fitting the high-z SNe 
distances. In §5 we derive the effects of lensing by a point- 
mass lens. This model is then generalised to a Hernquist 
profile in §6 and applied to observed galaxy-quasar lenses in 
§7. We present our conclusions in §8. 



2 THE TEVES EQUATIONS AND 

APPROXIMATIONS FOR GALAXIES 

Bekenstein's theory involves two metrics, one of which we 
denote by g M „ and is the metric in the Einstein frame, and 
the other is the physical metric which couples to matter, 
g^, following the notation of Bekenstein. The two metrics 
(and their inverses) are related by 

g^v = eT^g^ + (e~ 20 - e 20 )[/ M £/„, (3) 
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where <f> is a dynamical scalar field and C/ M = g^"U v is a 
dynamical, time-like vector field, normalised by g^UvU^ = 
1. Note that g M „ = g^ v in the limit that the scalar field 
<f> — 0. The dynamics of the vector field are governed by 
an action (cf. eq. 26 of Bekenstein), but for the problems 
involving lensing, which deal with either a static galaxy or 
a homogeneous universe, the vector field simplifies in the 
Einstein frame to 



= [(-<?«)" 5 c, 0,0,0], 



(4) 



i.e., a four-vector which is parallel to the time axis apart 
but with a different normalisation. The gravity sector of the 
theory is given by the Einstein-Hilbert action in the Einstein 
frame, 



5* 9 = J d 4 x - 



-99 



Pac g^ v 



(5) 



where Gk is a parameter of the theory, and here g is the 
determinant of the metric tensor, is the Ricci tensor in 
the Einstein frame, and Pac 2 is the energy density due to 
a cosmological constant. Matter is coupled to the physical 
metric by the usual action S m = fdx \/—gC m where C m 
is the Lagrangian for the luminous matter fields. Hence the 
physical energy momentum tensor 



[(pc 2 + p)u tl u„ + pg^ 



e*U u 



(6) 



where the renormalised physical four- velocity, u M , is treated 
as colinear with the time- like field, U^. The scalar field, cj> is 
governed by an additional action (cf. equation 25 of Beken- 
stein 2004). As a result, the scalar field <f> tracks the matter 
energy-momentum tensor distribution, satisfying an equa- 
tion, 



/' 



= 4nG K c- 4 [g^ + (1 + e-^)U»U v ] T„„, (7) 



where p(5^) is a function of 8^ , which is defined by 



(8) 



Here p and 5^ are non-dynamical fields, i.e., they are some 
fixed functions of the scalar field <j> an d the metric g^ v . As 
will become obvious ao can be identified with the accelera- 
tion scale in MOND. 

Near a galaxy the scalar field, 0, is quasi-static, hence 
we can neglect all time derivatives compared to spatial 
derivatives so that 



(0,V< 



(9) 



gtt 



Bekenstein argues that in the Einstein frame gtt 
-e 2 ' tN/c2 c 2 , with $iv(r) = $(r) - (f>(r)c 2 . Substituting 
these, and a simplified energy momentum tensor (equation 
6), into equation (7) for the scalar. Drop the pressure term, 
time derivatives and higher order terms, which are all small 
compared to the physical matter density term, p, the equa- 
tion now starts to resemble Poisson's equation, 



2V-V$^ 



2V- 



M rV6c 2 



8nGp, G = G K e' 2 ^,(W) 



where we setoff) « <f> c , which is the cosmological average of 4>, 
and define G as the physical gravitational constant. Hence 



$jv(r) is basically the Newtonian gravitational potential, 
and is related to the physical matter density by the New- 
tonian Poisson's equation. The equation of motion in these 
potentials is 



dt 2 



-V$(r) = -V<E>jv(r) 



(r)c 2 , 



(11) 



The above two equations link the gravitational potential $, 
the scalar field <j>, the matter density p and the motion dr/dt 
together. Plus a certain choice of the free interpolation func- 
tion, we fully specify the dynamics near a galaxy. As a first 
study we adopt a very simple choice of the free function with 



■£ 2 (i^) 2 
o 



< At < 1, 
otherwise . 



(12) 



As will become obvious, this choice of <5 2 is necessary to 
simplify the analytics of lensing by an extended galaxy. 

To suit our studies, here we have chosen a set of no- 
tations slightly different from Bekenstein. We avoid using 
Bekenstein's p, completely, which does not have the same 
role as Milgrom's p, which is a function of gravity g — \ V$|. 
To distinguish from Milgrom's p(g/ao) as well, we use /i(<5 2 ) 
to emphasise its dependence on the scalar field strength, not 
the overall gravity. * 

In our notations the toy model proposed by Bekenstein 
(his eq. 50) related <5j, and p by 



~2 



(i-A) 5 



(i 



(l-K/2)k jl \2 
8tt l-p.) 



(1 



(i-jf/2)fc p, > 

47T 1 — jl ' 



(13) 



where the parameters k and K are much less than unity. 
Fig. 1 compares these Bekenstein's toy model with our toy 
model. Our choice is different from Bekenstein's choice which 
allows a d 2 , with no upper limit and with a negative branch 
with subtle effects on cosmology. Nevertheless, both choices 
preserve the asymptotic relation p — > ^ in the weak grav- 
ity limit, which is essential for explaining galaxy rotation 
curves, the main success of MOND. We have also made lens- 
ing predictions for scalar field closely matching Bekenstein's 
toy model, this is given in the Discussion section. 



3 LENSING AND POTENTIAL IN TEVES 

3.1 Bending of Light Rays in Slightly Curved 
Space Time 

Light rays trace the null geodesies of the space time met- 
ric. Lensing, or the trajectories of light rays in general, are 
uniquely specified once the metric is given. In this sense 
light bending works exactly the same way in any relativistic 
theory as in GR. 

Near a quasi-static system like a galaxy, the physical 



* Our Gk is related to Bekenstein's G and K by Gk = G(l — 
K/2)- 1 , where K is a small dimcnsionlcss proportionality con- 
stant for the action of Bekenstein's vector field. Our p is re- 
lated to Bekenstein's scalar a by = 47rGxo" 2 . And 5^ = 
j~ 47r(l-K/2) l 2 . 



) 1 2 

-I y/3 where y is Bekenstein's y-parameter. The ao is 



related to the Bekenstein parameters k and I by ao 
S = (1 — /^/2)e- 2 ^. 



V3k 



and 
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Figure 1. shows the scalar field strength as function of ^ 
in Bekenstein's toy models (region shaded by dots with k and 
K between and 10). Also shown is our principal model (solid 
line cf. eq. 12), and the smooth model (dashed line at the lower 
boundary Bekenstein models, cf. eq. 80) 





Figure 2. Schematics of the lens geometry considered in this 
section. Top: Overall geometry. We indicate the angular diameter 
distance from observer to lens, D; and lens to source, Di s , the 
observed image angle 9, the bend angle a and radius of Newtonian 
bubble ro. Bottom: The geometry of a light beam past a point 
lens, indicating the position of a photon distance from the point- 
source, r, the distance of closest approach, b and the radius of the 
"Newtonian bubble", ro 



«()■ 



space-time is only slightly curved, and can be written in 
polar coordinates as 



2 , 2 

-c dr 



2*/c 



c 2 dP + e 



-2*/c 



df 



dl 2 = 



(14) 

(dr 2 + r 2 d<9 2 + r 2 sin 2 9dip 2 ). (15) 

where |$(a:)| < c 2 takes the meaning of a gravitational po- 
tential in a rectangular coordinate x = (xi,X2,x;i) cen- 
tred on the galaxy, we note that a non-relativistic mas- 
sive particle moving in this metric follows the geodesic 
w — V<1>(;e), which is the equation of motion in the 
non-relativistic limit. 

Consider lensing by the potential <l>(r) of a galaxy with 
a geometry as in Fig. 2. An observed light ray travels a 
proper distance l os = hs + loi from a source to the lens and 
then to an observer. Hence it arrives after a time interval 
(seen in polar coordinates by an observer at rest with respect 
to the lens) 



t ; 



-2S(r)/c 



2dl 



Jo c Jo 



2$(r) dl 



(16) 



where we used the fact that a light ray moving with a con- 
stant speed c inside follows the null geodesies dr — 0. As in 
an Einstein universe (cf. Bartelmann & Schneider 2001), the 
arrival time contains a geometric term and a Shapiro time 



delay term due to the $ potential of a galaxy. Gravitational 
time delay hence works as in GR, but with $ instead of the 
Newtonian gravitational potential, <l>jv(r). Note that we re- 
cover the GR-like factor of two in front of $. Images will be 
on extrema points of the light arrival surface, hence the fac- 
tor of two propagates to the deflection angle as well. Hence, 
to a good approximation, lensing by galaxies in TeVeS be- 
haves as in GR apart from different interpretations of the 
gravitational potential. 

For example, assume a spherical lens with a potential 
$(r). A light ray with an impact parameter b, we can reduce 
the eq. (109) of Bekenstein to 



a(b) 



dr 4b d<I>(r) 
^Vr 2 -& 2 dr : 



(17) 



where we take the weak field and thin lens approximation 
(i.e., we can drop higher order terms, assume the lens is far 
from the observer, the source li s oo and l oi — > oo, and 
approximating the closest approach r m i n ~ b). Interestingly 
this is twice the bending angle predicted from extrapolating 
non-relativistic dynamics (Qin et al. 1995). 

In fact, gravitational lensing in TeVeS recovers many 
familiar results of Einstein gravity even in non-spherical ge- 
ometries. For example, an observer at redshift z — sees a 
delay At ohB in the light arrival time due to a thin deflector 
at z — zi 



At bs(R-) 



D s 



2D,D, 



(R-R s ) 2 



J — i 



dl 



2$(R,Z) 



(18) 



as in GR for a weak-field thin lens, <I>/c 2 <§C 1. A light 
ray penetrates the lens with a nearly straight line seg- 
ment (within the thickness of the lens) with the 2-D coor- 
dinate, R — Di8, perpendicular to the sky, where Di(zi) = 
'oi/(l + zi) is the angular (diameter) distance of the lens at 
redshift zi, D s is the angular distances to the source, and 
Dis is the angular distance from the lens to the source. The 
usual lens equation can be obtained from the gradient of the 
arrival time surface with respect to R. 

Nevertheless, there are important differences between 
lensing in TeVeS and in GR. These are mainly in the pre- 
dicted metric for a given galaxy mass distribution, and the 
predicted metric and distance-redshift relation for the Hub- 
ble expansion, which we will come to. 

There are, however, conceptual differences between 
TeVeS and MOND. The modification of gravity in TeVeS 
is made through the factor V0| 2 c 4 aQ~ 2 ), hence depends 
on the scalar field gradient |V</>|, rather than |V$| as in 
MOND. Strictly speaking, and V-3? are generally not 
colinear except for special geometries, hence the two descrip- 
tions are not generally identical. 

In short Bekenstein's theory reduces to MOND in the 
non-relativistic and spherical limit. Later in the paper, we 
will work exclusively under the assumption that ~ 4> c ~ 
0, or e^ ~ 1. Hence we can drop the tilde sign without 
confusion. We will consider spherical systems only, where 
the magnitude of the gravitation field is given by equation 
(11), and ao is given by equation (1). 
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Figure 3. Upper panel: Angular diameter distance Da(z) (in 
units of the MOND critical distance Dq = c^a^ 1 ~ 6cii u ~ 1 ) 
vs. redshift , z, in two cosmologies: (Q m ,Cl\) = (0.04,0.46) 
("minimal- matter" cosmology, thick solid) and (0.25,0.75) 
(LCDM, thin dashed). The data are converted to angular diame- 
ter distances by multiplying by a factor (1 + z) and assuming 
Ho = 70 kms _1 Mpc~ 1 . We also label typical lenses at redshifts 
z = 0.1 - 1, and typical quasar sources at redshifts z = 1 - 3 (as 
labelled). Note Dj ^- Z ^ <g; 1 at all redshifts due to the low value 
of ao, which means that cosmological objects have too small an- 
gular diameter distances for any prominent deep-MOND effect in 
lensing, i.e., most of the bending is done in the Newtonian bubble 
of Fig. 2). Middle panel: The minimal- matter cosmology (TeVeS) 
fit to high-z SN la luminosity distance modulus is slightly poorer 
than the ACDM flat cosmology fit, although only by A\ 2 = 3.5 
or 1.9a; error bars indicate the ±lcr limit on the estimated Q\. 
Bottom panel: shows for models of different Q\ the sound hori- 
zon at 2 = 1000 (in degrees) in the above "minimal-matter" open 
cosmology and CDM-based closed cosmology. 



4 ANGULAR DIAMETER DISTANCES FOR 
LENSING IN A SIMPLISTIC 
MINIMAL-MATTER OPEN COSMOLOGY 

To describe gravitational lensing, it is essential to know how 
the background cosmology behaves in Bekenstein's TeVeS. 
Lensing requires knowing the luminosity distance Dl(z) and 
the angular diameter distance Da in a universe as functions 
of redshift. These can be predicted generally by 



(1 + z) 2 



sinh 



= D A (z) = 



Q 



1/2 



r 

Jo 



dz 
H(z) 



(1 + *) 



1/2 



(19) 



where Q.k is the usual curvature density parameter and H (z) 
is the Hubble expansion rate. To develop the cosmological 
model, we assume an isotropic and homogeneous physical 
metric in TeVeS (cf. equations 117 and 118 of Bekenstein); 

-c 2 dT 2 = -c 2 dP+a(t) 2 [d X 2 + sinhx 2 (d6 2 + sin<9 2 dV> 2 )] (20) 

for an open universe with the expansion factor a and the 
cosmological time t. 

We are interested in lensing in a reasonable cosmological 
model in TeVeS context, which should be in accord with 
Supernovae (SNIa) data. There are many degrees of freedom 
in TeVeS, which is beyond the exploration of this paper. One 
way is to consider a "mininal cosmology" in TeVeS context 
with just a low density baryonic universe plus some amount 
of cosmological constant model such that to fit SNe data. As 
argued by Bekenstein, the scalar field contributes negligibly 
to the Hubble expansion, with a ratio O(k) compared to 
the matter contribution, and k is small. Hence the minimal- 
matter cosmological model is largely the same as in the case 
of GR, and for latter discussion we set the energy density of 
the scalar field = (cf. Section VII, Bekenstein 2004). 
Under this crude approximation, we find that the Hubble 
parameter H(z) is given by 



H^(n m (i + zf + n A + nK(i + z) 2 ), (21) 



where Qk ~ 1 — — fl m - At low redshift the vacuum en- 
ergy is important and dominating. Choosing an appropriate 
vacuum energy, it is conceivable that one can make a TeVeS 
baryonic universe which expands at virtually the same rate 
as the CDM model. For example, we find a reasonable fit 
is an open cosmology with £7a ~ 0.46, Q m ~ 0.03, and 
Ho ~ 70km/s/Mpc. It is possible to fit the luminosity dis- 
tances of high-z SNe about as well as CDM up to a redshift 
of 1 - 2. Fig. 3 show the angular diameter distance as a 
function of redshift, for a best-fit minimal-matter cosmol- 
ogy and a standard ACDM model. Clearly both models are 
consistent with the high-z SNe distance moduli data set; 
the minimal-matter cosmology fit is slightly poorer than the 
ACDM fit, but only by 1.9a and is therefore admissible. 
This open cosmology has problems at very high redshift, 
e.g., the last scattering sound horizon will be too small. This 
undesirable feature might be an artifact of the our simple 
cosmology, rather than intrinsic to TeVeS; e.g., allowing for 
2eV neutrino Dark Matter to co-exist with baryonic matter 
in MOND would soften the problems (Sanders 2003). More 
detailed analysis is clearly needed. Despite its limitations 
our simplistic minimal-matter cosmology is sufficient for as- 
signing lensing and source distances in redshift z — 1 — 3. 
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Table 1. Lookup-table for our notations for lensing in TeVeS 
Parameters Meaning 



of eq. (18) with respect to R, and requiring images to form 
at extreme points, we find the lens equation 




: Oc/Oo 



Threshold for weak acceleration regime 
Gravity strength | ^7«J> | in TeVeS 
Newtonian gravity strength |V3>jv| 
Interpolation function 

Scalar field strength in quasi-static galaxies 
Angular diameter distance to lens 
Angular diameter distance to source 
Angular distance from lens to source 
Distance scale of weak acceleration 

Rescaled lens/source effective distance 
Newtonian bubble radius and angular radius 

Hernquist lens scale length and angular scale 

Angular size of closest approach to lens 

Deflection angle 

Reduced deflection angle 

Critical (Einstein ring) angular radius 

Rescaled angular radius of the lens 

Rescaled angular critical radius 



— — Dls ,, s 

a = -jj-a, a(b) = 



dr Ab d$(r) 
c 5 " Vr 2 - b 2 ' dr 



(24) 



Taking one more derivative to 9 S , we find that the amplifi 
cation A is given by 

d6 



A^ 1 = 



6d6 



(l-K-7)(l-K + 7). 



(25) 



5 MODELLING A SPHERICAL LENS AND A 
POINT LENS 

Having set out the basics of gravitational lensing and the 
cosmological model in the relativistic TeVeS theory, we now 
consider analytic solutions to simple spherical lensing mod- 
els. Equation (10) reduces to the MONDian equations for 
$(r), 



(22) 



if we assume the two gradients V$ and Vci are parallel 
(as in a sphere or far from a non-spherical body) or if the 
gradient is negligible (as in strong gravity). 

The simplest case of lensing is by a point-like mass dis- 
tribution. Chiu et al. (2005) have worked out the deflec- 
tion and time delay in the point lens case rigorously for 
TeVeS. We follow the more GR-like formulation of Beken- 
stein (2004), neglecting high order terms. The results of the 
two approaches are essentially the same when dealing with 
a galaxy potential. Qin et al. (1995) and Mortlock & Turner 
(2001) have previously worked out deflections and amplifi- 
cation by a point lens. Here we expand these earlier works 
by predicting divergence, shear and critical line for a general 
spherical lens. We will also fit the observed lenses later in 
the paper. A list of notations are given in Table I. 

Consider a geometry as illustrated by Figure 2, where 
a spherical lens at redshift zi bends the ray from a source at 
redshift z s . Following the convention in gravitational lensing 
to work with angles projected on the sky, we let the source 
be offset from the lens by an angle 9 S and form an image 
at angle 9, which is related to the physical length b (for the 
closest approach) by 



b_ 



(23) 



The spherical symmetry of the lens means that the line of 
sight to the lens, source and images lie in one plane. Many 
of the familiar results from gravitational lensing in Einstein 
gravity transfer directly over to TeVeS. Taking the derivative 



Here k and 7 are the convergence field and the shear field, 
still given by 

where K is the mean convergence within a circular radius. 

To offer more insight, we can re- write the deflection 
angle (cf. eq. 24) in terms of the photon position along an 
unperturbed path (the Born approximation) as 



j — 00 



2<?x(r), 



\A 2 +fc 2 , 



where I is the distance long the light path and 



9±(r) = 9(r)~, 



9(r) 



d<f>(r) 
dr 



(27) 



(28) 



is the gravitational force acting transversely to the direction 
of the photon motion at a distance r from the source. 

To simplify notations later on, it is helpful to scale the 
distances to a TeVeS distance scale Do defined by 

c 2 6c 

Do = — « ~ 25Gpc. (29) 
ao Ho 

Heuristically, this gives a characteristic distance in TeVeS 
at which a particle accelerated with a — ao would reach the 
speed of light (ignoring relativistic effects). Define r\ as a 
dimensionless number with 

V=y/4i,t=£,D=°!g!L. (30) 

L>o L>s 

The quantity r\ or £ characterises the geometry of the lens 
system, independently of the lens mass. For a point lens of 
mass M the meaning of r\ can also be recognised from the 

fact that rfi = = 0.8"^^=^^ is the 

conventional Einstein radius in the GR limit; this does not 
hold in TeVeS. Generally rj or £ represent the lens geometry 
and are independent of the lens mass M. 

With our choice of 8 function of jl relation (equa- 

tion 12) we can further simplify the gravity inside a spherical 
system. We find 



g(r) = V$ 



y/aog N 9n{t) < ao, 



(31) 



which is a special case of eq. (1) of Milgrom's theory. Here 
the Newtonian gravity 5jv(r) and potential $jv(r) are given 
by 

GM(r) _ d$iv(r) 



9n(t) = 



dr 



(32) 



where M(r) is the mass enclosed inside a spherical physi- 
cal (proper) radius r. In TeVeS, the gravitational field g for 
a spherically-symmetric source is made up of two contribu- 
tions, V0 from the scalar field and gjv from the Newtonian 
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gravity. For our choice of the free function we find that grav- 
ity in TeVeS is related to Newtonian gravity by 



V0c 



gN 



ro, 



Qn 



gN > ao, 
g N < ao- 



(33) 



5.1 Lensing by a Point mass 

Around a point-mass the gravity (cf. eq. 31) 



a (r /r) 2 , 
ao(ro/r) 



r < ro, 
r > r , 



(34) 



where we have defined a transition radius ro by j(r = 
ro) = ao, so that inside a Newtonian bubble, r < ro, 
around the point-lens we have strong Newtonian gravity 
where g oc 1/r 2 , and outside we have weak TeVeS gravity 
where g oc 1/r. At r = r we have 



GM GM 
ro = \ = — 5- = lOkpc 

a vl 



M 



v = (GMa ) 1/4 = 200 kms" 1 



10 n M© 
M 



10 n M, 



(35) 



(36) 



is the circular velocity in the weak gravity regime outside 
the Newtonian bubble. 

Substituting the TeVeS gravity g(r) (eq. 34) into 
eq. (27), we find the deflection angle 



4GM 
c 2 b 
2-kvI 



1_(A) +^sin- 1 A, o<r ,(37) 
V ro / c l ro 



b> r . (38) 



This result is straightforward to understand: for a light path 
with a large enough impact parameter b, we have r > ro 
everywhere along the line of sight so that the TeVeS gravity 
looks like that of an isothermal halo of circular velocity vo, 
and hence we recover the result of a constant deflection as in 
GR for isothermal halos. In the other limit where 6 ro, the 
line of sight will go through the strong gravity, Newtonian 
regime and the deflection approaches that of a point mass 
in GR, 



4t>g&o 4GM 



(39) 



c 2 b c 2 b 

This term also appears as the leading term for small dis- 
tances from the source in equation (38) . The extra terms in 
equation (38) are due to modifications as the light beam 
passes through the weaker, MOND-gravity regime, when 
r > ro- 

It is convenient to work with dimensionless quantities 
to find universal relations within TeVeS. Define 9o as the 
angular size of the Newtonian bubble with physical radius 
ro- 



ro 
Di 



GM 1 
ao Di 



= 2" 



M 



10 n M o 



We can then express the image angle 9 in terms of the di- 
mensionless angle y where 



do r 

and find the deflection angle satisfies 



(41) 




Figure 4. Shows dimensionless functions of the rescalcd impact 
parameter y = 8/80 in TeVeS (solid lines) and in GR (dashed 
lines): (a): deflection angle a/ (dor) 2 ); beyond the Newtonian bub- 
ble of angular radius 80, the deflection is a constant, (b): conver- 
gence k/?? 2 ; note the cusp at the transition between weak and 
strong MOND regimes; also k = in GR for a point lens, (c): 
shear 7/r; 2 ; the signal increases significantly within the New- 
tonian bubble. Panel (d) we shows two magnification profiles 
for a point-mass source in TeVeS (solid) and in GR (dashed), 
for values of the dimensionless MOND-lcns geometric parameter, 
v = 2^/1= 0.4 and = 1. 



a = J \ftr 
eon 2 



1 + sin y, 



y < 1, 
y> 1, 



(42) 



where rj is given in eq. 30. We can also transform the con- 
vergence and shear into dimensionless quantities. 



if I JL 



y < 1, 
y> 1, 




(43) 



(44) 



The amplification A — (1 — k— 7)(1— k +7)| _1 - These dimen- 
sionless results are plotted in Figures 4. We see that both 
the deflection angle and the shear decrease from the cen- 
tral point, just as for the point-like lens in Einstein gravity. 
But beyond the MOND-angle, 9 = 9o, the deflection angle 
is a constant, and the shear falls slowly as , just as it is 
for an isothermal sphere in Einstein gravity. This is perhaps 
not too surprising, since the aim of MOND was to mimic 
the rotation curve of a dark matter-dominated isothermal 
sphere with a point-like MOND source. As we have mod- 
ified gravity in the same way for gravitational lensing, we 
find a similar result. 

The shape of the convergence field, K./r} 2 , in Figure 4 can 
also be simply understood. For angles less than the MOND- 
angle, 0o, a light ray will pass for most of its path through 
the Universe in the weak MOND regime. We assume there 
is some large-scale cut-off to this and we return to the stan- 
dard cosmological model on very large scales. In this weak 
MOND-regime the light bundle will experience a conver- 
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gence. Then, when the light ray passes through the Newto- 
nian bubble, it experiences no further convergence, except 
for those rays that pass through the central point-mass. At 
the MOND-radius, 9 = 9o, there is a cusp as light rays inter- 
sect the edge of the Newtonian bubble. This cusp is a result 
of the discontinuity in g, given by equation (34), and other 
versions of TeVeS will give smoother transitions. For light 
rays that pass further away from the point-lens and never 
intersect the Newtonian bubble, the convergence falls off as 
y~ 2 , as for an isothermal sphere in Einstein gravity. The 
magnification pattern depends on the lens distance parame- 
ter r] or £, and is substantially different from GR prediction 
if r\ is not very small. 

Further insight into the behaviour of these results can 
be found by taking the limit y <C 1. In this limit we find 

1 



a 

K 

rj 2 

K 



1 2 

o9 



2 , 

V + 



4 4 

V V 




3 y 4 





2rf 
15 



y 



(45) 
(46) 
(47) 

(48) 



which show the first-order corrections to lensing by a point- 
lens in Einstein gravity. In physical coordinates the conver- 
gence field is 



2GMD U 



DuDi 



(49) 



DsDic* " w s ' s D S D 

which reduces to the usual result for a point lens in Einstein 
gravity when Do — > oo. Hence the first-order correction to 
the usual Einstein gravity results for gravitational lensing 

2 

for TeVeS is an added constant of 2£ = ^- . 

5.2 Critical line around a point lens 

In general, the critical radius is where the magnification di- 
verges. For a spherical lens, A — > oo when 9 S = 0. The crit- 
ical radius, 9 — 9 C , can be solved by setting 9 S = 9 — a = 
in the lens equation (cf. eq. 42). We find that the rescaled 
critical radius y c = |^ is a function £, satisfying 

»c<l, HO! 



77 



£ = — = 
s 4 



Vc 



4\/yl 

Vc_ 

2n 



' — 1 + 4 arcsin y c 
otherwise 



(51) 



which could be inverted approximately by interpolating the 
asymptotic relations, 



y 



-1/2 



(52) 



Note that we avoid the term "Einstein ra dius" , which nor- 
mally means an angle AGMDu/ D s Dic 2 in GR. The criti- 
cal radius 9 C = 9oy c (0 in TeVeS is also proportional to \[M, 
but depends on the distances through £ in a more sophisti- 
cated way. 

In the MOND limit for y c < 1 and £ < 1 we find the 
critical lines lie at 



(53) 



10.00 



1.00 r 



0.10 



0.01 




0.01 



0.10 



Lens Distance f 



Figure 5. Shows the critical radius y c = Di8 c /ro (rescaled by 
Newtonian bubble ro = \J GM/oq), as a function of the lens 
effective distance parameter, £ = ?? 2 /4 = DiDi a /D s /Do (rescaled 
by the MOND distance Do = c 2 /a ), where ao = 10~ 8 cms~ 2 . 
Crosses indicate 44 lens systems with data from the CASTLES 
lensing survey, including those in Table 2 and Tabic 3. The curves 
are for a point lens (solid), and a Hernquist profile lens with a 
scale length = ro/3 (dashed), and ro (lower dashed). 



In the limit that y c 3> 1, the TeVeS point-mass looks 
like an isothermal sphere in Einstein gravity, with a constant 
deflection angle. Here the magnification, convergence and 
shear fields are given by 



2k|" 



7 ; 



Any critical lines occur when k 



2tt£, or 9 C = 27t£0 o = 0.56" 



1/2 at a radius of 
~M D ls 



10 U M H D s 



(54) 



(55) 



In Figure 5 (solid line for point mass lens) we have 
plotted the dimensionless position of the critical lines, y c = 
9 c /9o, as a function of the lens geometric factor, £ = 77-^ • 
At small values of the critical line, and small values of £, we 
are in the strong-gravity regime and so we expect to recover 
the results of Einstein gravity, where y c oc In the weak- 
gravity MOND regime we expect a transition to equation 
(55), where y c oc £. Hence, we can use data from the observed 
measurement of critical lines around galaxies with known 
baryon content to test this prediction of relativistic-MOND 
theory: basically we compare the expected baryonic mass of 
a lens of observed luminosity L with the gravitational mass 
M inferred from the observed Einstein ring, which satisfies 
eq. (52). We carry out this procedure in §7.1. 



6 LENS MODEL FOR BARYONIC 
HERNQUIST PROFILE 

Early type galaxies are fairly well described by a nearly 
round distribution of light with a de Vaucouleur radial pro- 
file in projection. In this section we shall extend our results 
from a point-like lens to an extended, Hernquist-profile lens 
model. The latter has a density profile p(r) and enclosed 
mass profile M(r) given by 
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p(r) 



Mr h _ 

2-Kr(r + r h )' s ' 



M(r) 



r + r h 



-M, 



(56) 



where M is the total mass and is the core scale length. 
This Hernquist model is a good fit to the de Vaucouleur 
profile of elliptical galaxies, and is frequently used for quan- 
tifying these profiles (e.g. Kochanek et al 2000, Kochanek 
2003). Since in a TeVeS context the mass must follow light, 
a Hernquist profile should provide an accurate lens model. 





6.1 Gravity and scalar field 

The Hernquist model has a Newtonian gravity 
, , GM(r) GM 

fljv(r) = = ■ rj. 

r 1 (r + r h ) 2 

Hence in TeVeS the gravity is given by 
GM 



(57) 



9(r) 



(r + r h ) 
v 2 



r < r () - r h 



r > r - r h (58) 
r + r h 

where r — r — r h is now the radius of the Newtonian bub- 
ble between strong and weak gravity. Note that the gravity 
reaches a finite maximum at the centre of a Hernquist pro- 
file. We define 

e h = (59) 

as the angular size of the core length scale, r h , projected on 
the sky. A point mass model is obtained in case that 9 h — > 0. 

The scalar field gradient around a Hernquist profile 
galaxy is given by (cf. eqs. f2, 33), 

V "' :2 f°- r i r<n, -' V f , m 

otherwise. ^ ' 



ao 



= A(i 



f° 



+ rh 



1 — T ° 
r+r h 



Here the parameter jl = min(f , r ^ rh ) , hence satisfying 



jEi = min(f, —) = \± 



V0c 2 
ao 



(61) 



Integrating the above from r to the infinity, we find the 
scalar field 

0, 



(r) = 



%ln- 



GM 
c A (r+r h ) 



GM 



r < r - r h , 
otherwise. 



(62) 



Hence <f> contributes only beyond the Newtonian bubble. 

Figure 7 compares the gravity (solid curves) and Newto- 
nian field strength (plus signs) around a galaxy in a compact 
Hernquist mass model (th = 1 kpc) and an extended Hern- 
quist model (rn = 11 kpc). For our choice of ft, the scalar 
field contribution starts at small radii if the gravity g is be- 
low ao , it then peaks at a value ao /4 at a radius r = 2r — m , 
and then falls to zero again. 

Figure 8 compares the radial distribution of the circu- 
lar velocity for a point-mass and a Hernquist mass distribu- 
tion in non-relativistic MOND. In the point mass case we 
see that the velocity curve drops in a Keplerian fashion in 
the strong, Newtonian limit. At radii beyond the Newtonian 
bubble radius, the gravity enters the weak MOND-limit, and 
the velocity flattens to v = vo- The size of the Newtonian 
bubble scales with M 1 ^ 2 , and the terminal velocity scales 
with M 1 / 4 . The Hernquist models approach the point mass 



Figure 6. The geometry of a light beam passing a Hernquist lens, 
indicating the distance of a photon from the centre of the source, 
r, the distance of closest approach, b, the radius of the Newtonian 
bubble, ?*o, and the core radius of the Hernquist profile, r^. Top: 
illustrates situation when < ro- Bottom: illustrates situation 
when rj, > tq. 




Log r/kpe 

Figure 7. The rcscaled Newtonian field strength <?jv(r)/ao 
(crosses) and rescaled gravity g(r)/ao (solid curves) for an un- 
smooth function fi (cf. eq. 12) and a smooth function (cf. eq. 80) 
around Hernquist models of scale length rj, = 1 kpc or 11 kpc 
(thin or thick curves) for an elliptical galaxy M = 10 11 Mq. The 
scalar field contribution V<j>c 2 /ao is (g — gjv)/ao- 



model at large radii; but at small radii, the circular velocity 
curves are rising instead. 



6.2 Deflection by the Hernquist model 

We can now find the reduced deflection angle, a, for a light 
beam passing at a radius 6 — b/DL from the Hernquist lens 
centre. This can again be calculated from equations (27) and 
(28). Substituting in equation (58) we find 



a 

where 



OoOu + OoOhHo-u 
02 ~ e l 



(63) 



Hi 



where 



iin(j\/H) — arcsin(fc\/7f) 
!inh(j v / |if|)-arcsinh(fc v / | if l) 



H < 



(64) 



= 



otherwise. 



(65) 
(66) 



The different conditions in equations (64) to (66) correspond 
to different ray paths through the Hernquist lens, illustrated 
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.5e3 



■3 -2<i3 



.5e2 



1 

Log r/kpc 



Figure 8. Circular velocity for an elliptical galaxy (lower thick 
and thin lines) with M = 10 11 Mq for a point mass model 
(thin line) and a Hernquist model (thick line) with scale length 
= 3 kpc for an unsmooth function ft (cf. eq. 12) and a smooth 
function (cf. eq. 80); The boxes show the contribution of the scalar 
field in these two choices of fi. Likewise for a galaxy cluster with 
M = 10 13 Mq (crosses) for a point mass model and a Hernquist 
model with scale length = 30 kpc. The sharp break point is the 
transition radius ro from strong gravity to MOND weak gravity. 



by Figure 6. For instance, the upper ray on the top panel 
is outside the Hernquist scale length corresponding to \0\ > 
8h, ie H > 0, and is also outside the Newtonian bubble, ie 
\9\ > 0o — Oh- This latter condition also requires 0h < 0o, ie 
the Hernquist length is smaller than the Newtonian bubble. 

The lower ray on the top panel illustrates a different 
regime, where \0\ > 0o — 0h\ now the ray passes through the 
Newtonian bubble, corresponding to the condition given by 
equation (65). 

On the other hand, the bottom panel has 0u > 0o, ie 
the Hernquist length is larger than the Newtonian bubble. 
In this figure therefore, the condition in equation (66) always 
applies. Nevertheless, the two rays obey the conditions H > 
and H < respectively, as one ray passes through the 
Hernquist length and one does not. 

In summary, the conditions in equation (64) govern 
whether a ray passes through a Hernquist length of the lens; 
the conditions in equation (65) and equation (66) govern 
whether a ray passes within a Newtonian bubble which is 
larger than a Hernquist length. 

In the limit that Oh = 0, we find 



H 0u = H lu — — = — arcsinw 



A~ x = 



hd0 s 
0d0 



-(•-IK'-S) 



— a = (condition for Einstein ring) . In the latter case, 
the critical radius C = y c 0o is given by f- = 1. Substituting 
in eq. (63), we have 



J_ _ u c + H 0u y h y c 1 Hiy 
where yh = 0h/0o and 



= 



- Vh 
otherwise 



1 - Vh > \Vc 



(69) 



(70) 
(71) 



and Hjk is a function of H = 1 — ^ as in eq. (64). The 
above j/ c (£) relation can be used to weigh the lens mass M, 
and hence measure the M/L of the lens. Basically 



GM 
ao 



ro 



OcDi 

yc{i,yh) 



DiDi a ao 



(72) 



where the lens distance measure Di, £ and the critical line 
radius C are all observables. In the point mass limit we 
have y h — 0, and we find H = 1, u = \J max(0, 1 — jyl), 
Hiu = ■§ — arcsin(tt) = arcsin max(l, y c ), so the critical line 
satisfies a simpler relation (cf. eq. 50). Note the functional 
dependence of the rescaled critical radius y c and lens dis- 
tance measure £ is fairly complex even with our simplest 
choice of the free function /t (cf. eq. 12). It is therefore an- 
alytically challenging to work with more general functions 
for /i. 

Figure 5 illustrates how the critical radius increases with 
the distance parameter £ = if 2 1 A for various values of the 
Hernquist length scale ru; more concentrated models have 
bigger critical radii, with the biggest in the point mass limit. 
At small lens distances, we find that the critical radius of an 
extended Hernquist model can be very small; this is due to 
the fact that the Hernquist model has a surface density that 
diverges only logarithmically at small radius. For example, 
for a very extended Hernquist model with yh ^> 1 ^> y c , 

2 

we have u = H 0u = 0, H ~ — ^ < 0, rj~ 2 = H lu /y c ~ 



Vh 



1 arcsinh — . Hence 4 s - — > sinh 



-% < 0, 



v 

-1/2 



» i; 



10, |0|>6>o, 
which recovers the results for the point-mass (cf. eq. 45) . 

6.3 Amplification by a Hernquist model 

We can find the amplification A for the Hernquist lens from 
the lens equation 



the critical radius is relatively small for an extended Hern- 
quist model deep in weak gravity regime where (^yr^j <S 1> 

and Do = c 2 /a > A- 

The convergence k and shear 7 can also be derived an- 
alytically; the expressions are somewhat too lengthy to be 
produced here. Instead we illustrate k and the inverse of am- 
plification 1/A as a function of the impact parameter b for 
(•57) a hypothetical example of a Hernquist-profile galaxy clus- 
ter lens in Figure (9). Comparing with the point lens case, 
the convergence is much larger at small impact parameter 
for the Hernquist case. At very large radius, the Hernquist 
model approaches the point lens case. 



(68) 



where a = %fct is given by eq. (63). The amplification 
diverges when ^g- — (condition for radial arc) or S = 



7 COMPARING TEVES PREDICTIONS TO 
GALAXY LENS DATA 

7.1 Stellar mass of the CASTLES sample 

Now that we have developed models for lensing by point 
masses and Hernquist profile lenses, we will apply this for- 
malism to galaxy lens observations. In particular, we will 
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Table 2. Derived mass ratio M/M* for CASTLES double image lenses. Lens and source redshifts are listed, together with angular 
diameter distance to the lens and the rescaled lens geometry parameter £, and lens stellar mass (estimated from evolution and K- 
corrected I magnitude, except for two systems from R and K magnitudes as indicated by bracket). The lens mass M is calculated using 
our two techniques (image position and flux ratio) in TeVeS, for both the point mass model and the Hernquist model. For lenses composed 
of stars of normal stellar populations we expect M/M t ~ 1. 



Lens 

CASTLES 



Di 

(Mpc) 



M, 



(kpc) 



Ml, 



M 2 , 



M, 



M, 



Ml, Hern 



Q0142-100 

B0218+357 

HE0512-3329 

SDSS0903+5028 

RXJ0921+4529 

FBQ0951+2635 

BRI0952-0115 

Q0957+561 

LBQS1009-0252 

Q1017-207 

B1030+071 

HE1104-1805 

SDSS1155+6346 

SBS1520+530 

B1600+434 

PKS1830-211 

HE2149-2745 

SBS0909+523 



0.49 
0.68 
0.93 
0.39 
0.31 
0.24 
0.41 
0.36 
0.88 
0.78 
0.60 
0.73 
0.18 
0.72 
0.41 
0.89 
0.50 
0.83 



2.72 
0.96 
1.57 
3.61 
1.65 
1.24 
4.50 
1.41 
2.74 
2.55 
1.54 
2.32 
2.89 
1.86 
1.59 
2.51 
2.03 
1.38 



1190 

1400 

1560 

1010 

880 

730 

1070 

980 

1540 

1470 

1320 

1440 

540 

1430 

1070 

1540 

1200 

1506 



0.034 
0.013 
0.019 
0.032 
0.027 
0.022 
0.035 
0.027 
0.032 
0.032 
0.027 
0.032 
0.020 
0.028 
0.029 
0.031 
0.032 
0.019 



40.8 x 10 1() 

26.7 x 10 1() 
556. x 10 1() 
38.0 x 10 1() (-R) 
3.36 x 10 1() 
3.09 x 10 1() 
2.68 x 10 1() 
84.4 x 10 1() 
9.00 x 10 1() 
7.40 x 10 1() 
16.6 x 10 1() 
33.2 x 10 1() 
4.40 x 10 10 (K) 
28.0 x 10 1() 
4.00 x 10 1() 

14.8 x 10 1() 
20.0 x 10 1() 
655. x 10 1() 



1.6 



0.32 
0.29 
5.23 
0.80 
1.19 
1.50 
2.48 

1.32 



11.4 



0.34 

0.08 

0.0025 

0.7 

50.7 

0.75 

1.3 

1.1 

1.45 

0.53 

0.45 

2.1 

.43 

0.57 

1.9 

0.5 

0.5 

0.03 



0.19 
0.065 

0.53 

41. 

0.35 

1.38 

0.32 

1.04 

0.42 

0.23 

1.23 

0.33 
1.13 
0.94 
0.29 
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Figure 9. shows two models for a cluster of mass 1 X 10 14 M Q at 
D = 0.5Gpc with a Hernquist scale = lOOkpc (solid line) or 
= (dashed line). Panel (a): Lens convergence k. Panel (b): for 



(1 - k) 2 - 7 2 ] . The Hernquist 



le weak regime and can create 



the inverse of the amplification 
cluster model is only partly in th 
3 images very close to centre. A radial arc and a tangential arc 
can form near the critical points at b = 15kpc and 35kpc for the 
Hernquist model. Note that the arcs move inwards when going 
from a point model to Hernquist model. 



examine the consistency of the strong lensing predictions for 
galaxy lenses in the CASTLES survey (CfA- Arizona Space 
Telescope Lens Survey; Munoz et al 1999). 

Here we use galaxies with known double/quadruple 
lensed quasars, with critical lines estimated from the quasar 
separation, baryonic mass inferred from luminosity, and 
known redshifts from the CASTLES sample (Kochanek et 
al, 2000). In the spirit of TeVeS all of the lensing galaxies 
have M/L of order unity. However, to be rigourous, we must 
include the K-correction, the luminosity evolution with red- 
shift, and the possibility of significant gas and extinction 
from dust. 

Kochanek et al (2000, Table 6) have measured combined 
K-correction and evolution corrections for CASTLES lenses, 
and find a correction for the / band which varies with red- 



shift (0.36 < 2 < 0.88 as appropriate for most of our lenses, 
a low density open universe with flo = 0.3) corresponding to 
a ~ 11 percent mean offset in luminosity, which will not im- 
pact significantly on our conclusions about whether TeVeS 
provides reasonable mass-to- light ratios. 

To make the interpretation of our lens model more di- 
rect, we first estimate the stellar content of the high-z lens 
using its observed I-magnitude and a model for the spectral 
energy distribution of an old non-evolving stellar popula- 
tion. We first estimate the luminosity without K-correction, 
or evolution/reddening correction by 



D A (z)(l + z) 
10 pc 



10 



0.4(m 



.) 



(73) 



where Da(z) is the angular diameter distance to the lens, 
(1 + z) 4 gives the dimming effect, m A q is the Ao-band 
magnitude of the Sun. We then estimate the stellar mass by 
the simple formula 



M, 



: 7 



A 



xL, 7(1^) =7(Ao)10°' 4 ^,(74) 



where 7 (^7) is the estimated mass-to-light ratio of a typ- 
ical nearby elliptical galaxy observed at the emitting wave- 
length e.g., observations at I-band Ao = 8140.4 cor- 
responds to light emitting at z = 0.5 with a rest-frame 
5000A, approximately the V-band. The 



1+0.5 



wavelength 

parameter 71 parametrizes the countering effects of the K- 
correction, and the passive evolution/extinction after star- 
bursts; these effects tend to cancel each other for most cos- 
mologies and reasonable formation redshift of ellipticals, and 
models with very late/very early formation ellipticals give 
71 > and 71 < respectively (see Fig. 8 of Kochanek et al. 
2000). For the I-band we set 



4 1 71 
1 



7(8140.4) =4^, 

Lq 



71 « 0. (75) 
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This is calibrated using Fig. 32 of Worthey 1994, assuming 
a 12-Gyr old solar metallicity stellar population; while the 
I-band luminosity of a 12-Gyr old red elliptical galaxy is 
more than its V-band luminosity, it is actually comparable 
to the rest-frame V-band luminosity of its 5-Gyr old younger 
counterpart at z = 0.5. For three of our lenses without I 
magnitude, we use the nearest available band R, H, V or 
K magnitude. We make similar conversions of mass to light 
using Worthey models for these bands, and checked against 
the predicted redshift dependence of the lens galaxy colours 
V — I, R — I, I — H with those in the literature (cf. Fig.4 of 
Keeton et al. 1998, Fig. 1 of Kochanek et al. 2001). t 

Galaxies tend to be more dusty and gas rich at higher 
redshift. Falco et al. (1999) found that a rest-frame differ- 
ential reddening AE(B - V) ~ 0.07 m - 0.1 m on average 
for 23 lenses (mostly early-type galaxies). McGough et al 
(2005) have compiled extinction for 25 CASTLES lenses 
(see their Table 1); the extinction observed at I-band of 
most lenses is AAi ~ 8140 5 ^°+ z) -i x 3.1AE(B - V) ~ 
2(1 + z)AE(B -V) < 0.8 m , i.e., less than a factor of two 
correction of the luminosity. Note that we exclude five lenses 
with AE(B — V) > 0.4/(1 + z), which are mainly high red- 
shift late-type lenses. The consensus seems to be that gas 
and dust play insignificant roles in elliptical galaxies in gen- 
eral, and there is strong evidence against a large dust lane 
for a lens at z = 0.84 (Kochanek ct al. 2000b). In any case, 
only the far side of a spherical lens galaxy is affected by the 
dust lane, hence the dust effect on the total luminosity is 
likely mild, less than a factor of two. 



7.2 Critical lines and lens geometry test 

Now we compare the predicted position of critical radii in 
TeVeS with CASTLES data. Figure 5 shows the predicted 
scaling of the critical line opening angle, y c = 9 c /9o, with 
the lens geometry factor, £ = rj 2 /A. It is interesting that the 
range of galaxies in this survey all lie in the strong-gravity 
regime, where £ <C 1. This is due to the fact that the angu- 
lar diameter distance in a TeVeS universe has a maximum 
value which is much smaller than the Do scale of TeVeS (cf. 
Fig. 3). Therefore all of the lenses in the sample have a small 
£, and are never in the purely weak gravity regime. This 
" coincidence" is actually a natural consequence of the small 
value for do- has also been anticipated by Sanders (1999), 
who noted in passing that MOND effects become important 
below a threshold in the surface density of a galaxy cluster, 
while the minimal surface density of image-splitting galaxy 
clusters is always above this MOND threshold. Hence the 
deep-MOND effect is argued to be small. 

Secondly we note that the mean position of galaxies 
in the plot lie on the point-lens prediction. This appears 
to be a success for TeVeS, where the stellar mass can pro- 
duce reasonable size Einstein rings. On the other hand real 
galaxies are not point masses. Furthermore there are many 
galaxies which are outliers in this distribution; this scatter 





Figure 10. Upper panel: The TeVeS lens mass for Q0142 - 100 
using the image position technique; here we assume a point mass 
lens model. The two curves represent the permitted mass and 
source position for the two images; their intersection provides a 
unique consistent TeVeS mass for the object. Lower panel: The 
TeVeS lens mass for Q0142 — 100 using the image position tech- 
nique, with Hernquist scale length of 1.6kpc. 



/ 


- 


J 






Figure 11. The TeVeS lens mass for Q0142 - 100 using the flux 
ratio technique, with a point mass model. The curve represents 
the permitted mass and flux ratio for the two images; the straight 
line provides the measured flux ratio. The intersection provides a 
determination of the TeVeS mass for the object. Lower panel: The 
TeVeS lens mass for Q0142 — 100 using the flux ratio technique, 
with Hernquist scale length of 1.6kpc. 



in lens bending power is not expected in a TeVeS universe, 
as the displayed TeVeS predicted curve should be universal 
if galaxies can be approximated as point lenses. 

Thirdly the scatter of observed lenses cannot be ex- 
plained by extending the lens via a Hernquist model. As 
shown in Figure 5, Hernquist models have lower bending 
power, and so cannot explain the images with large separa- 
tion. 



7.3 Double image lenses 

We now conduct a more detailed study of double-image 
image lenses from the CASTLES survey in the context of 
TeVeS theory with the point and Hernquist models. Here we 
examine only lenses with almost co-linear double images, as 
non-co-linear images and quadruples cannot be modelled in 
detail by our spherical models. 



t We assume 7(20500/1) = 1 and 71 ~ 0.4 for K band, 

7(16000A) = 0.8 and 71 0.5 for H band, 7(5550A) = 7 and 

71 ~ 2 for V band, and 7(6500A) = 5 and 71 ~ for R band. The 
solar absolute magnitudes are from Allen's Astrophysical Quan- 
tities. These relations take account of dust in galaxies. 



7. 3.1 Source position method 

We determine the lens mass by two methods. The first is by 
inverse ray tracing, i.e. from one of the image positions we 
can predict the source position as a function of lens mass, 
using e.g. equation (24). By matching source positions from 
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both images, we therefore measure the mass according to 
the TeVeS/MOND theory. 

An illustration of this technique is given in Figure 10 for 
CASTLE lens Q0142- 100. This system includes two images 
at angular distances 1.9" and 0.4" from the lens, with source 
redshift 2.72 and lens redshift 0.49. In the figure we see that, 
for a TeVeS point-mass, we can calculate predictions for the 
real, pre-lensed angular position of the source as a function 
of mass for each of the two images; where these predictions 
intersect we have a unique measurement of the lens mass, 
which is here found to be M = (1.38 ± 0.002) x lO n M . 

The exercise can be repeated with a finite Hernquist 
scale-length for the lens. The bottom panel of Figure 10 
demonstrates that the resulting mass estimate is quite sen- 
sitive to Th\ for a scale length of 1.6kpc for our lens galaxy, 
the mass estimate is increased to (2.86 ± 0.004) x 10 n MQ. 
In order to find ru in each case, we use the half light ra- 
dius for the galaxy, r e as given for CASTLES lenses by 
Kochanek et al (2000) and using the relation of Hernquist 
(1991), Th — r e /1.82 to relate r e to rh- 



7.3.2 Flux ratio method 

The second method of measuring lens mass is to predict an 
amplification map on the image plane as a function of mass. 
This can be achieved by calculating the rate of change of 
deflection angle a with the image position 9 as a function of 
mass. From this one can calculate the rate of change of ap- 
parent position with respect to the unlensed source position 
9 S , hence the flux ratio of two images is predicted by 



Ai 
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1 



da i DlS 

ggl gs 

da 2 DlS 

ae 2 d s 



(76) 



Figure 11 shows this predicted ratio as a function of mass 
for the example double-image lens Q0142 — 100, for a point- 
mass TeVeS model; also shown is the empirical flux ratio 
between the two lenses. We find that the TeVeS mass con- 
sistent with the observed flux ratio is (7.8±0.09) x 1O 1O M . 
This is similar in magnitude to the mass found from the 
source position method above, but is clearly not consistent 
with it in detail. In order to obtain consistency, one requires 
a Hernquist model with r h — 2.3kpc. In contrast, the value 
of rn expected from half-light measurements is 1.6kpc, so 
modelling the extended nature of the lens mass cannot be 
used to account for all of the difference between the two es- 
timates. One may see this as a potential difficulty for TeVeS; 
alternatively we should note the simplicity of the model we 
are using (spherically symmetric, specific profile) which may 
account for some or all of the 43% difference between the two 
estimates. 

The bottom panel of Figure 11 shows the equivalent test 
with the measured r h — 1.6kpc; in this case we find a TeVeS 
mass of (2.52 ± 0.03) x 10 11 Mq. In this case, the mass is 
closer in magnitude to the TeVeS mass found via the image 
position technique (~ 13% difference), but is not consistent 
in detail; again, the simplicity of the lens model may be the 
cause of this. 



7.3.3 Results for CASTLES lenses 

We have applied the two TeVeS mass estimation techniques 
to a set of double image lenses from the CASTLES sur- 
vey (Table 2), for both point mass and Hernquist models. 
The lenses were selected to have known redshifts for source 
and lens, known magnitudes in F814, and to have only two 
images. We examine the TeVeS mass of these lenses in rela- 
tion to their absolute luminosity (from F814 magnitude) in 
Figure 12, for point lens and Hernquist lens cases with the 
position method. Seven of the lens galaxies do not have pub- 
lished half-light radii; for these we ascribed the median Hern- 
quist length of the galaxies with known radii, r h = 1.8kpc. 

We note that there is an apparent correlation between 
the TeVeS mass and the luminosity in the F814 band in 
both the point and Hernquist cases; this is broadly consis- 
tent with TeVeS, which predicts a strong correlation between 
mass and luminosity. However, there are some very signif- 
icant outliers, leading to a negative Pearson's correlation 
coefficient (-0.12 for point-mass; -0.04 for Hernquist). 

Finally we can derive the ratio of lens mass vs. stellar 
mass for each galaxy, using lens masses from either the posi- 
tion method or flux ratio method. In Figure 14 we show the 
mass-to-mass ratios M\ ens /M t from both of these methods, 
using the Hernquist model. 

Three points of interest arise from this plot. Firstly, 
we note that the mass-to-M» ratios calculated using the 
two independent methods closely agree, i.e. (M/M»)fi U x — 
(M/M t ) position- This shows that TeVeS is working well in 
terms of being a theory which appropriately includes gravi- 
tational lensing; if there were not a good level of agreement, 
it would suggest that different lensing predictions were not 
consistent in TeVeS. As it is, TeVeS is obtaining consistent 
results from lensing distortions of position and magnifica- 
tions. Nevertheless, for several lenses, e.g., SDSS1155+6346, 
the flux ratio cannot be reproduced with any lens mass, 
hence the empty entries in Table 2. This might reflect the 
fact that the flux ratios have been perturbed by microlens- 
ing or substructures. Overall the source position method is 
more reliable than the flux ratio method. 

The second point concerns the TeVeS mass and absolute 
luminosity for both point and Hernquist models. There is a 
strong correlation between the point and Hernquist mass es- 
timations themselves (see Figure 13). We can see that point 
lens models underpredict the required baryonic mass for a 
given set of image positions in relation to the Hernquist 
model, by a mean factor of 3.6. This confirms the value of 
extending our analysis to the Hernquist model. 

The third point of interest from this plot concerns the 
range of mass-to-M* ratios measured. All but two of the 
lenses are found to have M/M* between 0.5 and 2; this is a 
reasonably concentrated distribution. However, we note the 
existence of extreme outliers such as RXJ0921+4529 with 
M/L = 239 and HE0512-3329 with M/L = 0.16 (cf. Table 
2 for the mass ratio M/M* ) . 



8 DISCUSSION 
8.1 RXJ0921+4529 

As noted in the last section, there are several lenses with 
anomalous M/L. In particular the galaxy RXJ0921 with 
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Figure 12. Upper panel: TeVeS lens mass for double image CAS- 
TLES lenses using the image position technique (point mass), as a 
function of (converted) F814 absolute luminosity for these lenses. 
The lines indicate M/L = 4 (dashed). Lower panel: TeVeS lens 
mass for double image CASTLES lenses using the image position 
technique (Hcrnquist), as a function of F814 absolute luminosity 
for these lenses. Open circles show lenses where the median Hern- 
quist length has been used; for easy comparison, these lenses are 
also displayed with open circles in the upper panel, where Hern- 
quist length is irrelevant. Note five outliers indicated by five open 
diamonds. 
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Figure 13. Comparison of TeVeS masses obtained using the im- 
age position technique, for point masses and Hcrnquist models. 
The dotted line indicates equality of mass; open circles show 
lenses where the median Hcrnquist length has been used. 
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Figure 14. The lens mass to stellar mass ratio M/M* from two 
methods for Hcrnquist profiles using realistic length scales. The 
dotted line indicates equal mass-to-light ratios; open circles show 
lenses where the median Hernquist length has been used. The 
data are expected to be tightly clustered around M/M, ~ 1 in 
MOND/TeVcS. 



luminosity L = 0.84 x 1O 1O L0 requires either a point mass 
of 1.7x 1O 12 M or a Hernquist profile mass of 2.0 x f0 12 M Q . 
This lens is therefore required to have a M/L ~ 240, which 
could seriously challenge the TeVeS theory. 

Nevertheless, several explanations are initially available 
in defence of TeVeS for this system. The first possibility is 
that this lens is severely dimmed by obscuration of dust; 
however, the I — H colour of the lens appears normal, which 
would not be expected in the presence of dust. 

Murioz et al. (2001) also note an extended emission 
source B' offset from the B image for this system, and there 
is no extended counterpart in the image A. An important 
although unlikely explanation is that A and B are a pair 
of independent quasars with small separations and the same 
redshift, and the quasar B is off-centred from its host galaxy 
B' as well. A more plausible explanation is that B' is a mem- 
ber galaxy in the same cluster, projected by chance near one 
of the two split images of a lensed quasar. 

Another suspicion would be that this outlying mass is 
due to the use of too large a Hernquist length (as this is a 
lens where the median ru is used); however, the fact that 
the large mass is found even in the point-mass case allays 
this concern. 

Yet another possibility for the high M/L is the fact 
that the lens sits in the middle of an X-ray cluster; this 
will increase the lensing distortion associated with the sys- 
tem. This issue clearly requires follow-up using lens mod- 
els more sophisticated than those considered in this paper, 
which would be able to model multi-component lenses with 
different length-scales. 

Another system with a mildly high M/L is B1600+434; 
however, unlike our other lenses this is found to be an edge- 
on spiral galaxy. A Hernquist model will not provide a good 
fit to the baryonic matter in a spiral, so it is unsurprising 
that this M/L is an outlier. Indeed, with a point mass model, 
the mass-to- light is more reasonable (4.5 or 7.5). 

In summary, at least one outlier in our study succeeds 
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in overcoming some of the defences available to TeVeS; how- 
ever, the presence of a lensing cluster does call into question 
the level of disagreement between TeVeS predictions and 
data on this lens. 



8.2 Effects of different choices of jl and a 

One concern regarding the above result is that we used an 
unsmooth function jl (cf. eq. 61), which reduces the contri- 
bution of the scalar field to exactly nothing in the strong 
gravity regime. Since most lenses have impact parameters 
inside the Newtonian bubble, letting the scalar field con- 
tribute in the strong gravity could help towards the M/L 
discrepancy. Common choices of ft include 



'9/(9 + ao) 



9/V9 2 + a o 
. 1 - exp(-p/o ). 



(77) 



Another concern is that do could be slightly (by a factor of 
two in some systems) higher than 1.2 x 10~ 8 cms~ 2 given 
reasonable errors on disk galaxy rotation curve (RC) data 
(Sanders & McGaugh 2002). 

To address these concerns we consider a new smooth 
function jl so that 



V<j>c 2 
do 



(78) 



is a monotonic increasing function of jl, closely matching 
that chosen by Bekenstein (cf. eq. 13, as illustrated in Fig. 1). 
With a bit of algebra it can be shown that 



9 — 9n + V<^>c = g N + ^/g~N(Mj, 
and 
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(80) 



The gravity and rotation curve of a Hernquist galaxy or 
galaxy cluster in such a model are shown in Fig. 7 and 8. The 
curves are much smoother than those in our principal model. 
Here the scalar field contributes even inside the Newtonian 
bubble (strong gravity regime), unlike that in our principal 
model. Also this model and our earlier choice of jl bracket 
all common choices of jl (cf. Fig. 15). 

Interestingly, our new choice of jl also allows the deflec- 
tion of Hernquist lens model to be calculated analytically, 
and we find that the rescaled critical radius y c satisfies 
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This allows us to estimate the lens mass M by fitting the 
Einstein ring size; the solution has to be found by iteration 
because eq. 82 is an implicit non-linear function of M/L 
through y c and yn- 
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Figure 15. The MOND/TeVeS functions ft and 5$ (shifted down 
by a factor of 10). Common choices of ft are shown in symbols, 
which arc bracketed by our choices (thin and thick solid curves). 
All these choices converge asymptotically so that ft — + 8<p g/ao 
in the weak gravity regime (below the horizontal dashed line). 



Table 3. CASTLES lenses with anomalously large/small ^f^- 

Lens Comments 

RXJ0921+4529 2-image. Resides in cluster 

SDS1004+4122 2-image. Resides in cluster 

B1600+434 2-image. Edge-on spiral lens with dust lane 

MG2016+112 3-image. One or two lens planes 

B2045+265 4-image. Source near cusp caustic 

HE0512-3329 2-image. Gas-rich lens. 

SBS0909+532 2-image. Early- type, little gas/dust 

B0218+357 2-image. Lens is likely a spiral galaxy. 

RXJ1131-1231 4-image. V-magnitude. Cusp caustic. Elliptical. 

B1933+503 10-imagc components. Fits i? 1 / 4 -law. 

8.3 Constraints on jl from a larger sample 

Fig. 16 shows the application to a sample similar to those 
shown in Fig. 5, including all CASTLES double/four-imaged 
lenses with measured lens and source redshifts. Note the 
very large scatter overall. There are galaxies well above and 
well below the expected stellar mass, which cannot be sim- 
ply explained by extinction, evolution or K-correction. In- 
terestingly, SDSS1004+4122 is also in a galaxy cluster (from 
which the lens redshift was assigned) , similar to the situation 
of RXJ0921+4529. At the opposite end, the lens galaxy of 
HE0512-3329 is a gas-rich damped Lyman absorber galaxy, 
barely resolved at z = 0.93 with a size ~ 1 kpc compara- 
ble to the resolution of the Hubble Space Telescope (Gregg 
et al. 2000). Nevertheless, the stellar mass alone seems to 
exceed what is required for lensing. A similar situation per- 
tains to SBS0909+532, an early type lens without much gas 
and extinction (Lubin et al. 2000, McGough et al. 2005). 
Some features of these outliers are compiled in Table 3. 

Selecting only a well-behaved subsample (e.g., exclud- 
ing disky and/or dusty galaxies, galaxies with unknown half- 
light radius and/or unknown magnitude), we find that the 
scatter is still significant, with especially large deviations 
for three high redshift lenses. The Mi enB /M* shows mild de- 
pendence on the lens redshift zj, suggesting the need for 
K-correction. To make M lcns /M* ~ 1 for the z = 0.9 - 1.0 
lenses would require a model where 7(A Z ) increases by one 
decade (i.e., 71 = 2.5 m ) per unit redshift; this is too steep 
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even for purely K-corrected models where galaxies do not 
evolve (dashed line, taken from Fig. 8 of Kochanek et al. 
2000). While it has been suggested before that MG2016+112 
might involve two lenses at two redshifts (Nair & Garrett 
1997), a more recent model found that a single lens plane is 
sufficient (Koopmans et al. 2002). The early type (Sa) lens 
B2045+265 from the CLASS survey seems to be a robust 
case of a single lens with accurate NICMOS near IR data 
and radio data; the source is near a cusp caustic, and is 
split into three images on one side and one image 2" away 
on the other side of the resolved lens (Fassnacht et al. 1999). 
The peculiar flux ratios in both radio and near IR among 
the three images near the critical curve have been used to 
highlight the existence of dark substructures (Keeton et al. 
2003). 

Another way of describing the residue between mass 
and light is by using the theoretically uncertain modifica- 
tion function fi. We slide the function fi in between two 
expressions, say jli and fi2, which are that of equation (12) 
and that of equation (80) respectively by making the linear 
combination 
1 f 1 - f 

■= = ± + —r J -, -oo</<1. (83) 

The results with a slidable < / < 1 are shown as the 
upper thick part of the error bars in Fig. 16. Still we find 
M\ enB / M, =fc 1 for many systems. Model with / < seems 
necessary. 

Yet another possibility is that the value for do = a^ C = 
1.2 x 10~ 8 cms~ 2 from disk galaxy rotation curve (RC) mea- 
surements could be an underestimate. Larger values have 
been suggested in galaxy cluster studies (e.g., Pointecouteau 
& Silk 2005 and Sanders & McGaugh 2002). So suppose the 
theoretical value for 

a = (1, 2, 4, 8) x aa c = (1.2, 2.4, 4.8, 9.6) x 10~ 8 cms~ 2 .(84) 

Keeping ft as in eq. (80), adopting bigger ao has the effect 
of lowering M\ erla /M t (shown by the lower segments of the 
error bars). Nevertheless, a good match to Mi ens /M* ~ 1 
still seems difficult. A few systems seem to suggest an ao 
more than factor of 10 higher than the standard value ao 
from rotation curves of disk galaxies. In short, the recurring 
discrepancy between some lenses and MOND/TeVeS pre- 
dictions seems to go beyond the simplest lens models, and 
remains with many plausible functions of ft and plausible 
values for ao. 

A related point is that not all popular MOND p, corre- 
spond to a plausible Lagrangian for the scalar field. Zhao & 
Famaey (2005) noted that the standard function fl(g/ao) = 
s/\/ a o + 3 2 ma kes the gravity g (and all its functions) 
multi-valued for the same strength of the scalar field. This 
happens for all sharply increasing p,, including our princi- 
pal one (see Fig. 15). These models are perhaps undesirable 
since the scalar Lagrangian cannot be expressed uniquely by 
the scalar field strength. The lensing predictions of TeVeS 
are likely in a range narrower than depicted by the "error 
bars" in Fig. (16). 



9 CONCLUSION 

We have explored key properties of gravitational lensing 
as it occurs in Bekenstein's theory TeVeS, a fully covari- 
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Figure 16. Shows the lens redshift-dependence of the mass ra- 
tio Mi ens /M* derived from fitting the Einstein ring size of the 
CASTLES sample (circles). Circles with error bars are for well- 
behaved subsample of CASTLES lenses which have the Hernquist 
profile and little dust/gas. The upper thick part of an error bar 
brackets the range of uncertainty of the modification function ft 
(as shown in Fig. 15) for fixed ao = a^ c = 1.2 X 10~ s cms~ 2 . 
The lower part of an error bar with alternating thickness adopts 
a = (1 - 2)aJ c , (2 - 4)a^ c , (4 - 8)af c but fixes jj. as in cq. 80. 
Circles with central asterisks are lenses with unresolved scale- 
length (adopt 1.826^ = 0.1"). Circles with smaller circles are very 
faint lenses with unknown magnitude and scale-length (adopt 
/ = 21 m and 1.828 h = 0.1"). The solid horizontal line is the 
expected K-corrected value (unity) for plausible star formation 
and passive evolution of ellipticals. The dashed line shows the 
expected K-corrected value but assuming galaxies do not evolve. 



ant modified gravity theory. We have shown that TeVeS is 
frequently successful in predicting gravitational lensing phe- 
nomena; however, there are lens systems where TeVeS ap- 
pears to fail radically. 

We began by giving an account of TeVeS theory. Lens- 
ing is still governed in TeVeS by a gravitational potential $ 
which enters the metric in a similar way to General Relativ- 
ity. However, the interpretation of the gravitational poten- 
tial is different; rather than this being due to baryonic and 
dark matter, it is solely due to baryonic matter, which gener- 
ates a gravitational potential in excess of the GR prediction. 
The usual relations among lensing deflection angle, the lens 
equation, convergence, shear and magnification still hold ex- 
cept that their values as functions of the impact parameter 
have changed, as we have shown analytically for point-mass 
lens. E.g., the convergence of a point-mass lens is non-zero 
in TeVeS. Analytical lens models are also constructed for a 
Hernquist profile, which should be a good model for baryons 
in most elliptical galaxies. Conditions for generating arcs are 
discussed appertain to galaxy cluster scales. 

We then further tested the theory using galaxy lens data 
from the CASTLES survey. We find that the observed rela- 
tionship between critical lines and lens geometry (cf. Fig. 5) 
would be largely consistent with TeVeS if lenses were treated 
as point-lenses, but there are a handful of outliers. We also 
calculated TeVeS masses for CASTLES double-image lens 
galaxies, assuming point mass lenses or Hernquist lenses, 
and using image positions or fluxes to calculate mass. We 
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found that TeVeS/MOND provides an acceptable explana- 
tion for the lensing data in general. But in a handful of 
cases we obtain too large or too small mass-to-light ratios 
(cf. Fig. 14, Fig. 16, and Tables 2-3). Another way to put 
this is that we observe some cases of deviations from the 
expected baryonic content; or the ao derived from the criti- 
cal lines of some individual galaxies deviates markedly from 
the TeVeS/MOND expectation (cf. Fig. 16). More detailed 
photometric and spectroscopic data of these outlier lenses 
(see Table 3) are needed urgently to clarify the degree of 
the problem. 

On cosmology in the TeVeS context we noted the likely 
need for a cosmological constant (allowed in both TeVeS and 
GR). If Qa ~ 0.46 it is possible to construct a low density 
baryonic universe that agrees with current cosmological con- 
straints from supernovae and which equivalently provides an 
acceptable angular diameter distances at low redshifts. Such 
an open cosmology might not be the best for TeVeS since the 
sound horizon size at the last scattering is underpredicted. 
(cf. Fig. 3). 

This work has made initial tests of TeVeS using grav- 
itational lensing. Our analysis of TeVeS in the weak field 
regime suggests that our results might be valid for any rela- 
tivistic theory that asymptotes to MOND in the weak field 
regime. Further refinement of these tests (Shan et al. 2006) 
will involve examining non-circular TeVeS lenses to ascertain 
whether they fit the double and quadruple imaged lenses in 
more detail. Observed time delays between images of a dozen 
systems could also be constraining (Zhao & Qin 2005). Fur- 
ther work will also extend the analysis to a wider range 
of density profiles, including the beta-profile most suitable 
for analysing cluster lensing. In the meantime, we see that 
TeVeS succeeds in providing an alternative to General Rel- 
ativity in some lensing contexts; however, it faces signifi- 
cant challenges when confronted with particular galaxy lens 
systems. Finally, we see that gravitational lensing can be 
used as a useful approach to distinguish between theories 
of gravity, and to probe the functional form of the mod- 
ification function fi. In the mixed model (cf. Eq. 83) the 
gravity g rises more steeply with fl for models with smaller 
/ and even negative / (cf. Figure 15). E.g., the model with 
/ = (Bekenstein 2005) is more effective in bending light 
for the same baryonic mass than models with fi — x/(l + x), 
fl — x/vT+ x 2 , fl = min(l,a;) (cf. Figure 16). The model 
with / = 0.5 mimics the simple function fl — x/(l+x) in the 
intermediate regime. Any model with / > 0.5 is only only 
inefficient in lensing, but also always shows an undesirable 
peak in the scalar field strength 8^. This leads to unphysi- 
cal external field effect according to Zhao & Famaey (2005) . 
They also find that models with / < have a wide transition 
zone, which do not fit the sharply rising rotation curves seen 
in galaxies. Putting all these constraints together, it seems 
that the MOND/TeVeS models with < / < 0.5 have the 
best chance to survive. 
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